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ABSTRACT 

We present the results for the galaxy luminosity function as estimated from the 
Durham/UKST Galaxy Redshift Survey. This survey is magnitude limited to bj ~ 17, 
contains ^2500 galaxies sampled at a rate of one in three and surveys a ^-^4 x 
10^{h~^Mpc)^ volume of space. The maximum likelihood parameters for a stan- 
dard Schechter luminosity function are estimated to be M^ = —19.72 ± 0.09, 
a = -1.14±0.08 and 0* = (1.2±0.2) x 10"^ (/i^Mpc'^). Attempting to correct for the 
scatter in the observed magnitudes leads to a flatter faint end slope, a — — 1.04±0.08, 
which, combined with the different luminosity function shape, causes a higher nor- 
malisation to be estimated, (f>* = (1.7 ± 0.3) x 10~^ (/i^Mpc~^). Neither of these 
parametric functions provides a good formal fit to the non-parametric estimate of 
the luminosity function. A comparison with galaxy luminosity functions from other 
redshift surveys shows good agreement and the shape of the luminosity function now 
appears well-defined down to M^j ~ —17. There are some discrepancies between the 
different surveys for galaxies fainter than this absolute magnitude. However, our es- 
timate agrees well with that from the APM-Stromlo Galaxy Redshift Survey and we 
measure a fairly flat faint end slope. 



Key words: galaxies: luminosity function - galaxies: general 
tions - large-scale structure of Universe. 
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1 INTRODUCTION 

The galaxy luminosity function is a fundamental quantity 
in observational cosmology. Indeed, a knowledge of its form 
it is essential for an accurate analysis of galaxy clustering 
in magnitude limited redshift surveys (e.g. Efstathiou 1988) 
and also for a proper interpretation of the observed galaxy 
number-magnitude counts (e.g. Metcalfe et al. 1995a). On 
the theoretical side, the luminosity function is one of the 
first tests that models of galaxy formation and evolution 
must pass in order to be called successful (e.g. Cole et al. 
1994). 

Current interest in the observed galaxy luminosity func- 
tion takes the form of an accurate determination of the faint 
end slope and overall normalisation (e.g. Marzke et al. 1994; 
Vettolani et al. 1996), the relationship between the field 
and cluster luminosity functions (e.g. Driver et al. 1994; 
Kashikawa et al. 1995) and the evolution with redshift of 
the luminosity function (Lilly et al. 1995; Ellis et al. 1996). 
Meanwhile, theoreticians are attempting to improve the so- 
phistication of their modelling in order to predict more re- 
alistic luminosity functions (e.g. Frenk et al. 1996). 

The initial clustering results, redshift maps, etc. of the 



Durham/UKST Galaxy Redshift Survey were summar ized 
in the first paper of this series ( Ratcliffe et al. 1996| ). In 
this paper we present a detailed analysis of the luminosity 
function and space density of galaxies in this optically se- 
lected survey. We briefiy describe our survey in Section El 
The techniques used to estimate the luminosity function, ra- 
dial density and normalisation are outlined in Section H. The 
results from the Durham/UKST survey are then presented 
in Section kl Section M compares and discusses the current 
state of the galaxy luminosity function. Finally, in Section n 
we summarize our conclusions from this analysis. 



2 THE DURHAM/UKST GALAXY REDSHIFT 
SURVEY 

The Durham/UKST Galaxy Redshift Survey was con- 
structed usin g the FLAIR fibre optic system (Parker & 
Watson 199^ ) on the 1.2m UK Schmidt Telescope at Sid- 
ing Spring, Australia. This survey is based on the astrome- 
try and photometry from the Edinburgh/Durham Southern 
Galaxy Catalogue (EDSGC; Collins, Heydon-Dumbleton & 
MacGillivray 1988; Collins, Nichol & Lumsden 1992) and 
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was completed in 1995 after a 3-yr observing programme. 
The survey itself covers a ~20° x 75° area centered on 
the South Galactic Pole (60 UKST plates) and is sparse 
sampled at a rate of one in three of the galaxies to hj ~ 
17 mag. The resulting survey contains ~2500 redshifts, 
probes to a depth greater than 300/i~^Mpc, with a me- 
dian depth of ~150/i~^Mpc, and surveys a volume of space 
~4x 10''(/i-^Mpc)^ 

The survey is >75 per cent complete to the nominal 
magnitude limit of hj — 17.0 mag. This incompleteness 
was mainly caused by poor observing conditions, intrinsi- 
cally low throughput fibres and other various observational 
effects. In a comparison with ~150 published galaxy veloci- 
ties (Peterson et al. 1986; Fairall & Jones 1988; Metcalfe et 
al. 1989; da Costa et al. 1991) our measured redshifts had 
negligible offset and were accurate to ±150 kms~^. The scat- 
ter in the EDSGC magnitudes has bee n estimated at ±0.22 
mags ( Vletcalfe, Pong & Shanks 1995b) for a sample of ~100 
galaxies. This scatter has been confirmed by a preliminary 
analysis of a larger sample of high quality CCD photometry. 
All of these observational details are discussed further in a 
forthcoming data paper (Ratcliffe et al., in preparation). 



3 ESTIMATING THE LUMINOSITY 
FUNCTION, RADIAL DENSITY AND 
NORMALISATION 

We have estimated the galaxy luminosity function, 4>(.L), 
from the Durham/UKST Galaxy Redshift Survey firstly us- 
ing the very basic and common method of Schmidt (1968) 
and secondly using two maximum likelihood techniques; the 
parametric method proposed by Sandage, Tammann & Yahil 
(1979) and the non-parametric stepwise maximum likelihood 
method of Efstathiou, Ellis & Peterson (1988). These meth- 
ods have well defined error properties and the maximum 
likelihood techniques have been constructed such that they 
are unbiased by density inhomogeneities in the galaxy dis- 
tribution. They assume that the luminosity function has 
a universal form and hence the number density is separa- 
ble into a product of functions of luminosity and position, 
n{L, r) = (j>{L)p{r). Unfortunately, in using this assumption 
all the density and normalisation information is lost. To ob- 
tain the radial density information one uses a similar tech- 
nique to the above non-parametric method as proposed by 
Saunders et al. (1990). To determine the absolute normalisa- 
tion the minimum variance iterative technique of Loveday et 
al. (1992) is employed. This is a development of the method 
originally proposed by Davis & Huchra (1982). 

3.1 The Luminosity Function 

Probably the most basic estimator of the luminosity function 
is the 'l/Vmax' method due to Schmidt (1968) 



4,(L)dL = 



1 



ymax y-'-'i ) 



(1) 



where the sum extends over all galaxies in the L — » L ± dL 
luminosity interval and Vm,ax{Li) is the maximum volume 
that the galaxy of luminosity Li could be seen in (given the 
survey's physical and apparent magnitude limits). Unfortu- 
nately, this estimate is biased by density inhomogeneities in 



the galaxy distribution but is still useful for initial compar- 
isons. Error properties of this method can be estimated by 
a simple rms of the galaxies in each luminosity interval in 
question 



Yar{(l>) 



ymaxyJ-^i ) 



(2) 



For future reference this will be called the VMAX method. 
Methods that are independent of the galaxy density in- 
homogeneities can be constructed as follows. One can form a 
likelihood, jC, based on the probability of observing a galaxy 
of luminosity Li at redshift Zi in a magnitude limited red- 
shift survey 



Pi oc (t>{Li) / / (t>{L)dL 

where 



n 



Pi- 



(3) 



(4) 



The product extends over all of the A'' galaxies in the sur- 
vey and Lmin{zi) is the minimum absolute luminosity that 
a galaxy at redshift Zi could have and still be included in 
the survey. A maximum absolute luminosity could also be 
incorporated into equation |3|, altering the upper limit of the 
integral. In practice this makes little difference to any re- 
sults. The best estimate of ^(L) is then given when jC (or 
equivalently ln£) is maximised. 

In the parametric case one assumes a function form for 
(fiiL) and maximises equation W with respect to the pararn - 
eters of this assumed functional form (Sandage et al. 1979). 
I n keeping with tradition we use the 'Schechter function' 
( pchechter 1976 ) to describe 4>{L) 

0(L)dL = 0-(il)%xp(-i^)d(i_), (5) 

where the three parameters are a normalisation, 0*, a faint 
end slope, a, and a characteristic luminosity, L* , (or equiva- 
lently absolute magnitude, Af*). In practice 0* cancels from 
equation ^ and a maximum is determined in the {A4*,a) 
likelihood space. A further consideration is the effect of the 
observed scatter in the measured magnitudes. Efstathiou et 
al. (1988) modelled this by approximating these errors with 
a Gaussian distribution of zero mean and a rms. The ob- 
served (convolved) Schechter function, (j>o, then becomes 



MM) = 



2-n-a 



(j){M') exp 



-2^(^^' 



MY 



dM',{6) 



Given that equation H is a slightly naive approximation it is 
not obvious that it will substantially improve the quality of 
the fit. Finally, error properties of this method can be esti- 
mated by the deviations of £ from its maximum value. One 
can plot error ellipses in the (M* , a) plane by considering 
the contours corresponding to 



In/I = ln/I„ 



77X/3(n), 



(7) 



where Lmax is the maximum likelihood, n is the number 
of free parameters (namely two, a and M*) and /3 is the 
required confidence level for that number of free parameters. 
For future reference this will be called the STY method. 
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In the non-parametric case one assumes that </)(M) can 
be written as a series of constant steps across given luminos- 
ity intervals and one maximises to find the relative values of 
these steps ( pfstathiou et al. 1988 ) 

0(Af)=0fc, \M-Mk\<AM/2, k^l,...,Np, (8) 

where Np is the number of constant steps (or bins) of width 
AM. The maximisation condition (^f-^ = 0) quickly leads 
to an equation that can be solved iteratively 



EN 
1=1 



H(Mk-M^i„{zi))AM 



Y::."^^r 



■M^i„(zi))tl>jAM 



(9) 



where 



W{x)^ 


f 1 |a;| < AM/2 
1 otherwise ' 




and 






H{x) = i 


' X < 

^ 4- ^ ItI < 
2 + AM l-^l - 

. 1 X > 


-AM/2 
AM/2 
AM/2 



Finally, error properties of the 
the covariance matrix 



(10) 



(11) 



's can be estimated from 



Cov(<?!,fc) = 



d^\n£ 



(12) 



A constraint involving the (pk 's is usually introduced into the 
likelihood equation in order to ensure that the information 
matrix is non-singular and hence invertible as is assumed in 
equation [12. This constraint does not affect the shape of the 
maximum likelihood solution. We take a simple approach 
and fix one of the maximum likelihood (jtk's to be constant. 
We the n assume that one c an neglect the off-diagonal ele- 
ments (Saunders et al. 1990). Using this assumption gives 



Var((^fc) 




In practice the cross-derivatives are small and therefore our 
approximation will only slightly underestimate the errors in 
the (f)k- For future reference this will be called the SWML 
method. 

A major problem with the STY method is that it will 
always return a maximum likelihood solution regardless of 
the assumed parametric form and how well it represents the 
actual luminosity function. Therefore, it is necessary to test 
the goodness of fit of this solution with a likelihood ratio 
test. We assume that the non-parametric SWML method 
provides a good representation of the actual luminosity func- 
tion shape. Let £i be the likelihood calculated using the 
maximum likelihood solution of the given parametric form 
(STY) and let £,2 be the likelihood calculated using the max- 
imum likelihood solution of the (j>k's (SWML). Efstathiou et 
al. (1988) have shown that 21n(yCi/£2) approximately be- 
haves like a x^ statistic with {Np — 3) degrees of freedom. 
However, for the answer to be independent of the bin size 



and the number of bins, the likelihood £1 should be calcu- 
lated using the SWML likelihood formula (which produced 
£2) with the set of (/)fe's from 



J<j}{L)dN{L) J(l,{L)L^dL 



(13) 



JdN{L) jLldL 

The integrals in equation O are over the luminosity interval 



in question, [L^ — ^ 
Saunders et al. 1990). 



Lk + 



Ail 
2 



(Efstathiou et al. 1988; 



3.2 The Radial Density 



Using a similar method to Section 3.1 one can estimate the 
radial density function, p{r). First consider the probability 
of observing a galaxy of luminosity Li at redshift Zi in a 
magnitude limited redshift survey. As before a likelihood is 
formed from the product of these probabilities 



«=1 Jm 



p{zi) 



a\z 



;„(i.)] 



p(zi 



, dz ) 



dz 



(14) 



where Zmi-n(Li) and Zmaa;(iO are the minimum and max- 
imum redshifts at which a galaxy of luminosity Li could 
be seen and still be included in the survey. Saunders et al. 
(1990) assume that one can write pij) as a series of step 
functions in radial distance and then solve for the maxi- 
mum likelihood p^'s by iteration as before. Finally, error 
properties of the pfc's can be estimated from the appropri- 
ate covariance matrix and we use similar approximations to 
those made in Section 3.1 to simplify the analysis involved. 



3.3 The Normalisation 

By their method of construction these maximum likelihood 
techniques cannot provide an overall normalisation. There- 
fore, one must adopt other techniques to estimate the mean 
space density of galaxies, n, and the normalisation of the 
luminosity function, 0*. 

The expected distribution of the number of galaxies as 
a function of distance r is given by 



n{r) ^ fnV{r) S{r), 



(15) 



where / is the sampling rate of the survey, n is the mean 
spatial galaxy density of the survey and V{r) is the volume 
element of the survey at a distance r. The radial selection 
function of the survey, 5(r), is given by 



S{r) 



r 



41^1. 



^(r)] 



4>{L)dL 



c„ '^(^)'^-^ 



(16) 



where Liow is some minimum abs olut e galaxy luminosity 
and Lmin(r) was defined in Section 



3.1 



If a Schechter func- 
tion is assumed then the above integrals become incomplete 
Gamma functions, r(a-|-l,x). The mean spatial galaxy 
density is also related to the luminosity function by 



4){L)dL, 



(17) 



and so using the luminosity function from either equation H 
or o one finds a relationship involving n, (jj* and an integral 
over the shape of the assumed functional form. Using this 



4 A. Ratcliffe et al. 



formalism Davis & Huchra (1982) showed that an unbiased 
estimator for fi is given by 



£;^-^ S{r)w[r)dV' 



(18) 



where w(r) is any weight function we choose and \rmin,rmax\ 
is the distance range over which we use galaxies. One choice 
is to weight all galaxies equally, whereby to = 1. However, 
Davis & Huchra (1982) showed that the minimum variance 
in n occurs when 



■]{r) = - 



1 



where 



J3(r-c 



+ 4tv fnJ3{rc)S{r) 



X ^{x)dx, 



(19) 



(20) 



Tc is the scale on which J3 converges to a maximum value 
and ^{x) is the galaxy 2-point correlation function. Given 
that this weighting depends on the quantity we are inter- 
ested in (fi), Loveday et al. (1992 ) h ave developed an it- 
erative scheme involving equations hs and ha. This method 
should produce the minimum variance estimate of n if J3(rc) 
converges on a sc ale Vc smaller than the survey. The variance 
of n is given by ( Davis & Huchra 1982 ) 



Var(n) = -^ '-^ -^-1^ '- . (21) 

In practice the final answer for n (and hence (^*) depends 
little on the value of 3-3,(rc) used. Therefore, we adopt the 
value 47rJ3 (re) = 5000/i~'^Mpc'^ and doubling or halving this 
only makes a few per cent difference to n. Also, we are free 
to choose the value of Lio^ (or equivalently Miow) in equa- 
tions hq and O as long as wc are consistent and throw away 
all the fainter galaxies from the sum in the numerator of 
equation hq. Obviously n is very sensitive to the value of 
Miom used but <j)* is very stable over a range of ~5 magni- 
tudes. 



Table 1. STY maximum likelihood results and normalisations 
assuming the pure Schechter function of equation H, (f>, and the 
convolved Schechter function of equation H, ij}o ■ 







Mi 



a 



1) 



Prob. 
0* (/i3Mpc-3) 



-1.14 ±0.08 

-19.72 ±0.09 

0.16 

(6.4 ±1.1) X 10-2 

(1.2 ±0.2) X 10-2 



-1.04 ±0.08 

-19.68 ±0.10 

0.22 

(4.4 ±0.8) X 10-2 

(1.7 ±0.3) X 10-2 



where 

dL{z) = {l + z)r{z), (25) 

is the luminosity distance in h^^Mpc and kcorr is the k- 
correction. A simple k-correction is used for all galaxies in- 
dependent of morphological type 



f^corr — fClZ + K2Z 



(26) 



where fci = +3.15 and k2 = -0.29 ( [Broadbent 1994| ). Given 
the redshift range of interest this k-correction is more than 
adequate for this analysis. 

The survey's magnitude limits were chosen by the au- 
thors to maximize depth and minimize observational incom- 
pleteness. In this case each of the 60 UKST fields has a dif- 
ferent magnitude limit and our best sample contains 2055 
redshifts with {mum) ~ 16.86 ± 0.25 mags and an average 
completeness of 75 per cent. Selection effects like this have 
been sho wn to cause no systematic biases in our methods of 
analysis (Ratcliffe 199(:). One can check the completeness of 



this sample by using Schmidt's (1968) (V/Vmax) test. For 
a uniform distribution of galaxies in a complete magnitude 
limited sample (V/Vmax) = 0.5, with a rms dispersion of 
l/\/T2iV (in the absence of clustering). If (V/Vmax) is sig- 
nificantly lower or higher than 0.5 then we are missing ob- 
jects at high or low redshift, respectively. For the above best 
sample we find that (V/Vmax) = 0.50 ± 0.01 and so we are 
not systematically missing galaxies at any redshift. 



4 RESULTS FROM THE DURHAM/UKST 
GALAXY REDSHIFT SURVEY 

We use galaxies in the distance range [5,350]^-^Mpc, where 
the minimum distance comes from requiring a reliable red- 
shift distance estimate (relatively unaffected by peculiar ve- 
locities) and the maximum distance is due to the magnitude 
limits of the survey. We assume a go = ^, A = cosmology 
and so comoving distances are given by 



r{z) 



\HoJ 



VTT^ 



(22) 



where Hq = lOOh kms^^Mpc-^ is the Hubble constant and 
c is the velocity of light in kms-^. Volumes are calculated 
from the comoving distances and use 

dn 



V 



-r'{z), 



(23) 



where dfl is the solid steradian angle of the survey. Absolute 
(M) and apparent (m) magnitudes are related by 



5 Ig 0(1,(2) = m — M — 25 — kcorr{z), 



(24) 



4.1 The Luminosity Function Shape 

Using the magnitude limits from the previous section the 
STY parametric solution has been calculated from the 
Durham/UKST survey for a pure Schechter function (equa- 
tion H) and a convolved Schechter function (equation a). 
Following Metcalfe et al. (1995b) we use a = 0.22 in equa- 
tion M for the rms scatter in the observed magnitudes. Small 
variations in this parameter do not significantly affect the 
final results. The maximum likelihood results for a and Al* 
are shown in Table hi and assume h = 1. The SWML non- 
parametric solution has similarly been calculated from the 
Durham/UKST survey, taking ~20 iterations for 5 s.f. con- 
vergence. These maximum likelihood solutions are all shown 
in Fig. nl along with the basic VMAX estimate of the lumi- 
nosity function, where we have scaled the 0's to agree in the 
bin contaiiung the most galaxies (A/^j — —19.75). We are 
only considering the shape of the luminosity function here 
and so the absolute normalisation is still arbitrary. The in- 
set of Fig. hi shows the joint 68 per cent error ellipsoids for 
the two STY solutions as calculated from equation [71 while 
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Figure 1. The galaxy luminosity function maximum likelihood solutions for a stepwise function (dots), a pure Schechter function (solid 
curve) and a convolved Schechter function (dashed curve). Also shown for comparison is the basic VMAX luminosity function estimate 
(crosses). These four solutions are scaled to agree at Mj,^ = —19.75 but the overall normalisation is arbitrary at this stage. The inset 
shows the STY maximum likelihood results and joint 68 per cent error ellipsoids on both parameters in the (M*, a) plane. 



the errors quoted in Table hi are the la error on an individ- 
ual parameter. The VMAX and SWML error bars on Fig. gj 
were calculated from equations Q and y_3[ respectively. 

The observational incompleteness described in Section 
is not explicitly accounted for in the luminosity function 
shape analysis presented here. To get around this problem 
one can weight the Pi 's of equation S by the appropriate com- 
pleteness rate depending on the UKST field and apparent 
magnitude interval in question. In practice, this correction 
made very little or no difference to the either the STY or 
the SWML maximum likelihood solutions. 



A comparison of the two STY solutions in Fig. hi shows 
that the main effect of the magnitude errors on a luminosity 
function of this shape is to pull (;/>(Af) down at faint mag- 
nitudes while pushing it up slightly at bright magnitudes, 
essentially flattening it overall. The VMAX and SWML es- 
timates both agree well except at the very bright and very 
faint ends. Such a faint end excess is expected if there lo- 
cal overdensity, see Fig. H, and rejecting all galaxies within 
r < 20fe~'^Mpc brings the VMAX estimate into line with the 
SWML one. The STY and SWML estimates both have rel- 
atively flat faint end slopes with no convincing evidence for 
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conv 



100 200 

r (h-'Mpc) 

Figure 2. The maximuin likelihood estimate for the radial den- 
sity function. The solution has been normalised to unity in the 
observed distance range as indicated by the dotted line. Large 
fluctuations of order ~50 per cent are present on ~50/i~^Mpc 
scales. Also, the dashed line shows a simple straight line minimum 
X^ fit to the data points. This confirms the visual impression that 
the radial density is falling with distance. 




100 200 

r (h 'Mpc) 

Figure 3. The observed galaxy number-distance histogram, 
N(t), from the Durham/UKST survey using the magnitude limits 
described in Section W. The smooth curves show the expected dis- 
tribution for a homogeneous sample with this survey's magnitude 
limits and completeness rates using the maximum likelihood pure 
Schechter function (solid c urv e) and convolved Schechter function 
(dashed curve) of Section 4.1 



a faint end upturn down to Mb, ~ —15. The results of the 
likelihood ratio test described in Section 



3.1 



are also given 

in Table ry. It is seen that neither of these Schechter func- 
tions provide good formal fits to the SWML function, with 
the convolved Schechter function being marginally better. 
However, both agree well with the general SWML shape on 
a qualitative level. 



4.2 The Radial Density Shape 

The maximum likelihood estimate of the radial density func- 
tion has been calculated using the method outlined in Sec- 

The error bars on p(r) were 



tion 3.2 and is shown in Fig, 
also calculated using the techniques described in Section 3.2 
The non-parametric stepwise solution converged to 5 s.f. 
after ~20 iterations. This solution has been normalised to 
unity over the plotted distance range. 

The observational incompleteness previously described 
is not explicitly corrected for in this method. One can again 
weight the pi's implicit in equation [14 by the appropriate 
completeness rate to account for this effect. In practice, this 
correction makes only a small difference to the estimated 
p{r). 

Fig. shows that fluctuations in the observed radial 
density function are of order ~40-60 per cent and occur on 
~50/i~^Mpc scales. As we will see the radial distances of the 
'peak and trough' fluctuations agree well with the observed 
number-distance histogram. The observed large local over- 
density for r < 20/i~^Mpc probably causes the slight faint 
end excess seen in the VMAX estimate of the luminosity 
function. Also, it is interesting to see that a straight line 



minimum x^ fit to the radial density function (dashed line) 
is significantly different from the expected flat slope (dotted 
line). (This difference remains significant even if we remove 
the first point from the fit.) Specifically, the < 200/i~^Mpc 
region is ~20 per cent overdense and the > 200/i~^Mpc re- 
gion ~20 per cent underdense, the error being ~4 per cent 
on each measurement. This is discussed further in Section |5[ 

Fig. M compares the observed galaxy number-distance 
histogram with two homogeneous models. The magnitude 
limits and completeness rates from the previous section are 
used for models and data alike. The two models are the 
STY maximum likelihood solutions from Table hi for a pure 
(unconvolved) Schechter function and a convolved Schechter 
function. They have been normalised to match the total 
galaxy number in the sample. These homogeneous models 
do not provide a detailed fit to the data because of the ef- 
fects of strong galaxy clustering on scales > 20/i~'^Mpc. It 
can be also seen that both models overpredict the numbers 
of galaxies at r > 200/!.~^Mpc, with the convolved Schechter 
function model being slightly worse in this respect. 

Comparison of the observed Nir) from Fig. H with the 
radial density function from Fig. H shows that the 'peak 
and trough' fluctuations in p{r) agree well with the 'spikes' 
in N{r). As a consistency test, in Fig. we plot the prod- 
uct of the two homogeneous N{r) models with the observed 
p{r) and compare with the observed N{r). The agreement 
is impressive, implying that the initial assumption of the 
luminosity function having a uni vers al form (allowing a sep- 
aration of variables, see Section 3.1) is justified. 
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Figure 4. The observed galaxy number-distance histogram, 
Af(r). from the Durham/UKST survey (shaded histogram, see 
Fig bl). The two bold histograms are the product of the 
Durham/UKST radial density function (see Fig. p|) and the two 
number-distance model curves (see Fig 



4.3 The Normalisation 

The results of the iterative estimator of equations y_8| and |l9| 
for n are given in Table hi, with ~5 iterations needed for 5 
s.f. convergence. The value of 0* is estimated from numeri- 
cal evaluation of the integral in equation h? Use of a simple 
u; = 1 estimate for n (instead of equation 19) makes a ~5 
per cent difference to n. The formal error from equation Ell 
(using the weighting in equation M) gives a 7 per cent er- 
ror in n. This can be compared with a 10 per cent error for 
w — 1. However, the dominant error comes from the uncer- 
tainty in the selection function (due to the uncertainty in 
the Schechter parameters M* and a) and is ~15 per cent. 
This gives a total error in the normalisation of ~17 per cent. 
Also, for this analysis we choose Miow = — 15 and as previ- 
ously mentioned in Section 3.3 varying this parameter only 
causes a small (~ few per cent) change in (j)* . 

Obviously for estimates of n and (ji* the effects of obser- 
vational incompleteness must be corrected. Therefore, in this 
analysis the sum in the numerator of equation uM is weighted 
by the appropriate completeness rate of the UKST field and 
apparent magnitude interval of the galaxy in question. 

One can comment on the ~40 per cent change in <j!>* 
observed with the two different parametric forms. This is 
caused by a combination of the different relative luminosity 
function shapes (and the corresponding integrals over them 
in equations la, u/^ and uM and the different Schechter func- 
tion parameters (particularly a). This is discussed further in 
Section ^ 
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Figure 5. Comparison of the normalised Durham/UKST lu- 
minosity function with those estimated from recent galaxy red- 
shift surveys. All the magnitudes have been transformed to the 
blue bj system. This entails shifting the Las Campanas Gunn- 
r magnitudes by the average galaxy colour, (fej — r)^ = +1.1 
(Tucker 1994; Lin et al. 1996), while the CfA Zwicky magnitudes 
are shifted by bj - Mz = -0.45 (Shanks et al. 1984). All of 
the luminosity functions plotted here are corrected for magni- 
tude errors by convolution of a Gaussian error function with a 
Schechter function. For reference, we also plot the pure (uncon- 
volved) Durham/UKST Schechter function. 



5.1 Comparison with Luminosity Functions from 
Other Galaxy Redshift Surveys 

Table ti shows a comparison between the parameters of 
some recently completed galaxy redshift surveys (Loveday 
et al. 1992; Marzke et al. 1994; Lin et al. 1996) and the 
Durham/UKST survey. Note that the comparisons given 
here are for convolved Schechter functions because they are 
the preferred fits as quoted in the literature. However, as 
was seen in Section 4.1, while the convolved Schechter func- 



tion did provide a slightly better fit to the SWML lumi- 
nosity function, the significance was marginal given the ex- 
tra complexity added to the analysis. As was mentioned in 
Section 4.2 the pure (unconvolved) Schechter function pro- 



5 DISCUSSION 



vides a more realistic fit to the observed number-distance 
histogram. These normalised luminosity functions are plot- 
ted in Fig. ra where magnitudes have been transformed to 
the blue bj system, albeit with only a naive magnitude off- 
set. The Las Campanas galaxies are offset by the average 
galaxy colour, namely {bj — r)^ = -1-1.1 (Tucker 1994; Lin 
et al. 1996). The CfA galaxies are offset by bj - Mz = 
—0.45, which comes from the comparison of galaxy number- 
magnitude count data in Shanks et al. (1984). 

Looking at Table Q shows that the Las Campanas and 
CfA surveys contain considerably more galaxies than the 
Durham/UKST or APM-Stromlo surveys. However, in spite 
of this, all of these surveys sample similar sized volumes 
(within an order of magnitude) due to the different observ- 
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Table 2. Comparison of some recent galaxy redshift survey parameters and their corresponding maximum likelihood (convolved) 
Schechter luminosity functions. While the original passband of the survey is quoted in the left-hand column, the J\/*'s in the right-hand 
column have been transformed to the bj passband so that a direct comparison can be made {{hj — r)^ = -f 1.1 and bj — Mz = —0.45, 
Tucker 1994; Lin et al. 1996; Shanks et al. 1984). 



Survey 


Passband 


^Sal 


Durham/UKST 


bj 


-2100 


APM-Stromlo 


bj 


~1800 


Las Campanas 


Gunn-r 


~ 19000 


CfA 


Zwicky 


~9000 



Volume (h ^Mpc '•^) 

4 X 10^* 
1 X 10^ 
1 X 10'' 
1 X 10"^ 



-1.04 ±0.08 
-0.97 ±0.15 
-0.70 ±0.05 
-1.0 ±0.2 






0* (?t3Mpc-3) 



-19.68 ± 0.08 (1.7 ± 0.3) x 10^^ 

-19.50 ±0.13 (1.4 ± 0.2) X 10-2 

-19.19 ±0.02 (1.9 ± 0.1) X 10-2 

-19.25 ±0.3 (4.0 ± 1.0) X 10-2 



ing strategies used. The accuracy of the quoted <^*'s reflects 
both the number of galaxies and, more importantly, the total 
volume surveyed. 

Comparing the results in Table H and Fig. n shows that 
the general features of these luminosity functions agree well 
down to Mb J ~ —17. Namely, that the characteristic abso- 
lute magnitude at the 'knee' is M^ — —19.4, the faint end 
slope is flat, a ~ —1.0, and the overall normalisations agree 
well (apart from the CfA one which may be biased high by 
local inhomogeneities) . However, fainter than Mtj ~ — 17 
one does find some discrepancies in the values of the faint 
end slopes. The Durham/UKST and APM-Stromlo surveys 
agree well and retain a reasonably flat slope. Marzke et al. 
(1994) claim that the CfA survey has a signiflcant excess of 
galax;ies in this region, a ~ —1.3. Lin et al. (1996) measure a 
declining faint end slope, a ~ —0.7, from the Las Campanas 
survey. 

We can comment on the results from these other sur- 
veys. Firstly, the CfA excess comes mainly from a very local 
region (r < 25h~^Mpc), which is probably not a represen- 
tative sample of the Universe. Therefore, while an excess 
might exist, its significance could be in doubt. Secondly, 
not only is the Las Campanas survey selected in a differ- 
ent passband (which would cause different relative fractions 
of galaxy types to be observed), but also this survey is bi- 
ased against observing low surface brightness galaxies (be- 
cause of observational selection effects). Therefore, if there 
is any correlation between low surface brightness and in- 
trinsically faint absolute magnitudes then one can explain 
their deficit. We also comment on the preliminary luminosity 
function results from the ESO Slice Project (ESP; Vettolani 
et al. 1996). They measure a normalisation approximately 
twice that of the above surveys and a rising faint end slope, 
Q ~ —1.2, for galaxies fainter than Mi,j ~ —17. At present 
the significance of their faint end excess is unknown and so 
we cannot comment on that result. However, given that this 
survey is complete to a fainter magnitude limit than those 
discussed here, bj ~ 19.5, one might expect a higher nor- 
malisation because of the steep slope of the observed galaxy 
number-magnitude counts in this range (e.g. Metcalfe et al. 
1995a). 

We also note that the form of the Durham/UKST and 
APM-Stromlo luminosity functions are in good agreement 
with the combined colour-dependent galaxy luminosity func- 
tion estimates of Shanks (1990) and Metcalfe et al. (1997). 
Thus, the increasing slope of the luminosity function of blue, 
late-type galaxies cancels with the decreasing luminosity 
function slope of the redder galaxies, leaving an overall lumi- 



1.5 - 



0,5 




100 200 

r (h 'Mpc) 

Figure 6. Comparison of the radial density functions, p(r), esti- 
mated from the Durham/UKST and the APM-Stromlo surveys. 
These two surveys map a similar overlapping Southern region of 
space, with the Durham/UKST survey sampling a ~1500 sq. deg. 
area and the APM-Stromlo survey a ~4300 sq. deg. area. These 
two solutions have been normalised to unity within each survey. 
The dashed line shows the minimum x^ straight line fit to the 
Durham/UKST data from Fig. |. 



nosity function which has a flat faint end, which is observed 
here. 

Finally, it is possible to be too pessimistic about the 
statistical robustness of luminosity function estimates from 
local magnitude limited surveys (e.g. Driver & Phillipps 
1996). We emphasise again that the Durham/UKST lumi- 
nosity function presented here agrees very well with that 
of the APM-Stromlo survey, while the observed differences 
in the other surveys could be caused by slight systematic 
problems. 



5.2 The Radial Density Distribution 

Fig. |6| gives a comparison of the radial density functions, 
p{r), for the Durham/UKST and APM-Stromlo surveys. 
These surveys overlap each other in a Southern region 
of the sky but map the space slightly differently. The 
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Durham/UKST survey maps a 'wedge' (~20° x 75°) with 1 
in 3 sparse sampling, while the APM-Stromlo survey maps 
a thicker 'wedge' (~55° x 80°) with 1 in 20 sparse sampling. 
These two solutions have been normalised to unity within 
each survey. It can be seen that there is good agreement be- 
tween the two observed p(r)'s despite the fact that they do 
sample slightly different regions of space. In particular, the 
'peak and trough' fluctuations are observed in roughly the 
same places for both surveys. Also, the falling radial den- 
sity seen in the Durham/UKST survey again appears in the 
APM-Stromlo survey. 

It has previously been suggested that a very large local 
void, > 100^~^Mpc in extent, could explain the low normal- 
isation and steep slope seen in the bright galaxy number- 
magnitude counts {bj < 17), without the need for evolution 
(e.g. Shanks 1990). The observed trend in p(r) lends support 
to the idea that large-scale structure can affect the counts 
over a wide magnitude range. However, it should be noted 
that the decreasing trend in p(r) seen in Fig. n extends over 
the whole range of these surveys, out to 2: ~ 0.1. Since the 
counts beyond bj = 17 rise to have a 2x higher (p* than at 
brighter magnitudes this means that p{r) will have to rise 
sharply just beyond the 2 ~ 0.1 range of the present sur- 
vey if a large scale structure explanation of the steep galaxy 
counts is to be possible. 

The alternative possibility that the steep counts slope 
is caused by galaxy luminosity evolution is constrained by 
the fact that, at bj < 21, the galaxy M* seems unevolved, 
which confines any evolutionary explanation to affecting 
only the intrinsically faint galaxies (Broadhurst, Ellis & 
Shanks 1988). The observed decrease of p(r) with distance 
means that any such evolutionary increase in the numbers of 
faint galaxies has to be even larger in the range < 2: < 0.1 
than previously expected to explain the increase in galaxy 
counts at 17 < fej < 19. 

We have also considered the possibility that the de- 
creasing radial density p{r) could be caused by the effect 
of cosmology on the volume element. As noted by Weinberg 
(1972), the effect of qo on the volume element with redshift 
is first order, dV oc (1 — ^qoz)dz; thus even sd, z — 0.1 
there is a 15% change in volume for a unit change in go- 
However, the p{r) displayed in Fig. n was calculated using 
go ~ 0.5. Assuming a lower value of go such as go = 0.05 or 
go = —0.5 would therefore increase the volume element at 
high redshift and decrease the density at high redshift still 
further, acting in the opposite manner to that required to 
reduce the redshift dependence of p{r). Thus only a signif- 
icantly higher value of go than go = 0.5 would produce a 
more homogeneous result and we regard this as less likely 
than the alternative large-scale structure or evolutionary ex- 
planations. 

5.3 Comparison with the Cluster Galaxy 
Luminosity Function 

All of the galaxy luminosity functions mentioned so far have 
been estimated from the field. We can qualitatively compare 
our results (and those from other redshift surveys) with es- 
timates of the galaxy luminosity function from clusters. 

A faint end excess in the cluster galaxy luminosity func- 
tion (a ~ —1.4) has been known to exist in the Coma and 
Virgo clusters for quite some time (e.g. Abell 1977; Metcalfe 
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Figure 7. Comparison of the normalised Durham/UKST pure 
Schechter luminosity function with the cluster galaxy luminosity 
function model of Driver et al. (1994) and the galaxy formation 
model of Frenk et al. (1996). The model of Driver et al. (1994) 
has been made to agree with the Durham/UKST results at M* 
and 4>* . For the Frenk ct al. (1996) results the normalisations are 
fixed by the model. 



1983; Binggeli, Sandage & Tammann 1985). Recent work has 
confirmed these initial results and extended the observations 
to other clusters (e.g. Biviano et al. 1995; Kashikawa et al. 
1995). Driver et al. (1994) and Wilson et al. (1996) have also 
studied clusters at higher redshifts (z ~ 0.2) and observe a 
similar steeply rising upturn at the faint end. The overall 
conclusion is that a single Schechter function cannot provide 
an adequate fit to such luminosity functions. In all of this 
recent work the cluster galaxy luminosity function appears 
better fit by a 2-component Schechter-like function with a 
flat bright slope (q ~ — 1) and a steep upturn (a ^ —2) for 
galaxies fainter than Al* + 2. 

In Fig. W we compare t he S TY results for the pure 
Schechter function of Section 4.1 with the model proposed 



by Driver et al. (1994). This model approximates the gen- 
eral cluster luminosity function discussed above and has 
been translated in the Mt j and (j) direction to match the 
Durham/UKST results in Af * and 0* . As might be expected 
(given the above discussion) this form does not give a par- 
ticularly good representation of the Durham/UKST results 
at the faint end of the luminosity function (similarly for the 
APM-Stromlo and Las Campanas luminosity functions). On 
the other hand, it is tempting to say that the results from the 
CfA and ESP surveys match this form of luminosity func- 
tion. Overall, the errors in the faint ends of both the clus- 
ter and fleld luminosity functions are probably large enough 
such that no significant inconsistencies currently exist be- 
tween the two functions. 
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5.4 Comparison with Current Models of Galaxy 
Formation 

One can also qualitatively compare our results with those 
from current models of galaxy formation and evolution. 
These models approach the galaxy formation problem by 
combining the results of N-body/hydrodynamic simulations 
with semianalytic modelling techniques (e.g. Cole et al. 
1994). 

In order to successfully model galaxy formation one 
needs to include many distinct physical processes, as well 
as an underlying cosmological structure formation model. 
These processes include the evolution of dark matter halos, 
the dynamics of gas cooling, star formation and feedback, 
the evolution of the stellar populations that form and how 
galaxies merge and interact. One approximates these pro- 
cesses by a set of simple rules which then form the basis of a 
semianalytic model of galaxy formation (e.g. Kauffmann et 
al. 1993; Lacey et al. 1993; Cole et al. 1994). The models are 
then used to predict and compare with the observable prop- 
erties of the galaxy distribution, such as the faint galaxy 
number-magnitude counts, galaxy colours and the galaxy 
luminosity function. 

Fig. M also compares the pure Schechter function STY 
results from Section Wjj with the semianalytic galaxy forma- 
tion model of Frenk et al. (1996). The details of this fiducial 
model are described in Cole et al. (1994) with recent updates 
from Frenk et al. (1996) and Baugh et al. (1996). One can see 
that, given the complex modelling involved, this is a reason- 
able approximation to the luminosity function around M*. 
Note that the excess of very bright galaxies has been previ- 
ously documented (e.g. Frenk et al. 1996). Also, the model 
predicts a steep faint end slope to the luminosity function. 
Although this is more in line with the faint end excess seen 
in the CfA and ESP surveys, we have argued that these are 
less reliable estimates of the local galaxy luminosity func- 
tion, which we believe is better represented by the results 
from the Durham/UKST, APM-Stromlo and Las Campanas 
surveys. 



6 CONCLUSIONS 

We have estimated the galaxy luminosity function from the 
Durham/UKST Galaxy Redshift Survey. We use standard 
maximum likelihood techniques that are unbiased by den- 
sity inhomogeneities and an optimal weighting function for 
estimating the normalisation. For a standard Schechter func- 
tion we find M^j = -19.72 ± 0.09, a = -1.14 ± 0.08 and 
(j)* = (1.2 ± 0.2) X IQ-^h^Mpc-^. 

Correcting for the observed scatter in our measured 
magnitudes causes a fiatter faint end slope to be mea- 
sured {a = —1.04 ± 0.08), although this does not signifi- 
cantly improve the quality of the fit. A combination of the 
change in the shape of the functional form and the fiat- 
ter a then causes a higher normalisation to be estimated, 
(j)* = (1.7 ±0.3) X 10"^/i^Mpc-^. 

Comparison with other recent estimates of the galaxy 
luminosity function from redshift surveys gives good agree- 
ment for absolute magnitudes brighter than Mtj ~ —17. 
Fainter than this absolute magnitude there are some dis- 
crepancies in the measured value of a, ranging from -0.7 to 



-1.3. The Durham/UKST survey lies approximately in the 
middle of this range and is in very good agreement with the 
results from the APM-Stromlo survey. 

The radial density function has been estimated from the 
Durham/UKST survey using a maximum likelihood tech- 
nique. This shows evidence for a falling galaxy density with 
radial distance, suggesting that large-scale structure may af- 
fect the form of the galaxy number magnitude counts in the 
range 14 < 6j < 19. 

We have shown that the result of forming the product of 
the radial density function with homogeneous models which 
assume the above luminosity functions gives good agreement 
with the observed number-distance histogram. This result is 
consistent with the initial assumption that the galaxy lumi- 
nosity function has a universal form. 

Our luminosity function results for the field agree with 
those derived from galaxies in clusters, except at the faint 
end where the cluster luminosity function appears to be 
steeper. 

Finally, we have compared our results with the predic- 
tions of current models of semi-analytic galaxy formation 
and evolution. These also disagree with our luminosity func- 
tion at the faint end, where the galaxy formation models 
predict a significantly steeper luminosity function slope. 
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